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Abstract 

Lead clusters deposited on Si(lll) substrates have been studied at low temperatures using scanning tunneling microscopy 
and spectroscopy. The current- voltage characteristics exhibit current peaks that are irregularly spaced and varied in height. 
The statistics of the distribution of peak heights and spacings are in agreement with random matrix theory for several clusters. 
The distributions have also been studied as a function of cluster shape. TJjiX . 
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I. INTRODUCTION 

In the nearly 10 years that have passed since the first 
observation of discrete energy levels in metallic clusters 
there still remains the question of how the levels are sta- 
tistically distributed in these systems. It has been sug- 
gested that random matrix theory (RMT)i is applica- 
ble to the statistical properties of the spectra of metallic 
clusters in much the same way that it is applicable to 
the slow neutron resonant spectra^ observed in the 1950s 
and 1960s. However experimental verification of the ap- 
plicability of RMT in these systems is still challenging 
because of the difficulty in gathering a sufhcient number 
of levels to analyze their statistical distribution. 

The first observations of discrete energy levels or "par- 
ticle in a box energy levels" in metallic clusters were made 
by Ralph, Black and Tinkham^ in 1995. The energy lev- 
els were observed as irregular steps contained within the 
Coulomb staircase in the current-voltage characteristics 
of clusters that were fabricated using a fixed tunneling 
geometry with metallic electrodes. These uneven steps 
in the current-voltage measurements might be a conse- 
quence of random matrix theory (RMT).>A This sugges- 
tion arises from earlier predictions that address these 
systems from various theoretical standpoints, including 
Efetov's supersymmetry derivation.^ While these expec- 
tations are theoretically well established they are exper- 
imentally difficult to realize because of non-equilibrium 
effects^ and capacitive charging energy terms that have 
a tendency to mask the energy levels in mesoscopic sys- 
tems. Subsequent tunneling experiments have been per- 
formed on metallic clusters^ and semiconducting dotsS'— 
which have yielded results similar to those of Ralph et 
al. , although the nature of the level statistics still remains 
elusive. 

In contrast to metallic clusters, experimental work has 
exhaustively addressed the issue of distributions of level 
spacings and eigenfunctions of quantum dots fabricated 
from two dimensional electron gas systems of various pre- 
defined shapes.— In these systems, electron-electron in- 



teractions dominate transport and the level spacing dis- 
tributions appear to be Gaussian, while the distributions 
of the amplitudes of the eigenfunctions follow a Porter- 
Thomas distribution which is a signature of random ma- 
trix theory (RMT). 

Random matrix theory (RMT) and quantum chaos 
were merged in the conjecture put forward by Bohigas, 
Giannoni and Schmitii in 1984. This conjecture states 
that the nearest neighbor energy level spacings of clas- 
sically chaotic systems should be distributed according 
the Gaussian Orthogonal Ensemble (GOE), or Wigner- 
Dyson^^ distribution and this conjecture is strongly sup- 
ported by aggregated numerical studies. 

The Wigner-Dyson distribution which describes the 
statistical distribution of nearest neighbor energy levels 
normalized to the mean energy level has several impor- 
tant features as described by Porterii^ First, the proba- 
bility of having nearest neighbors with zero spacing dis- 
appears. Second, the probability of a level spacing is 
linear in energy before approaching a maximum, with 
the maximum occurring close to the mean energy and 
the tail of the distribution being fairly small. This is 
in contrast with completely random levels (a classically 
non-chaotic system) where the distribution is Poissonian. 
For the latter distribution the probability is largest at 
zero level spacing. The absence of small spacings in the 
Wigner-Dyson distribution is known as the: "repulsion of 
energy levels" This is the key ingredient of a Wigner- 
Dyson distribution and distinguishes classically chaotic 
from classically non-chaotic systems. 

Besides the distributions of the eigenvalues which are 
most often addressed experimentally, it is also possible 
to study the statistics of the amplitudes of the eigen- 
fimctions within the context of RMT. The corresponding 
distribution is called the Porter-Thomas^'^ distribution. 
The Porter-Thomas distribution is simply a statement 
about the amplitude of a wavefunction at any given 
point is a random variable. The distribution of the square 
of a random variable, , which is Gaussian distributed, 
is the Porter-Thomas distribution. 
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It is also interesting to note that it is possible, in some 
systems, to tune a transition from regular to chaotic be- 
havior in which the distributions of the nearest neighbor 
level spacing change from Poisson to Wigner-Dyson. This 
has been found numerically in the study of the Rydberg 
levels of the hydrogen atorn^ in which the transition is 
tuned by the application of an increasingly strong mag- 
netic field and experimentally in diamagnetic helium by 
tuning the excitation energy in the presence of a mag- 
netic field.— 

Here we present experimental results that indicate 
that statistical distributions of highly irregularly shaped 
Pb clusters follow RMT. This is an extension of earlier 
worki^ and its purpose is to address in detail experimen- 
tal issues and present additional results which led to our 
conclusions. The present paper is divided into five major 
sections as follows: in Section II we treat the device ge- 
ometry that was configured to study the discrete energy 
level spectra of metallic clusters. In Section III, we re- 
late energy spectra to the geometry and discuss in detail 
the observed features of the spectra. In Section IV we 
discuss the results of the preceding section in terms of 
statistical distributions of both the eigenvalues and the 
amplitude of the moduli of the eigenfunctions. In Sec- 
tion V, we demonstrate the use of scanning tunneling 
spectroscopy to resolve real time images of a quantity 
that is proportional to the square of the amplitude of the 
eigenfunctions in these systems. 



II. DEVICE GEOMETRY 

Before describing the current- voltage characteristics, it 
is worthwhile to summarize some of the salient features 
of the device configuration relevant to this work that 
are different from other spectroscopy measurements on 
metallic clusters. The tunneling geometry that was used 
is different from that used previously in three important 
respects. First, Pb clusters were grown on a semicon- 
ducting substrate. This substrate was highly resistive 
with a resistivity at room temperature greater than 1000 
ri-cm.ll Second, the clusters were fabricated using the 
buffer layer assisted growth technique that was developed 
by Weaver and co-workersi^ in which Xe, an inert gas, is 
used as a buffer layer to control the size distribution of 
the clusters. This physical technique allows the clusters 
to land softly onto the substrate and involves no use of 
chemicals or organics in the growth process. Third, fab- 
rication of the clusters and spectroscopic measurements 
were carried out in situ in an ultrahigh vacuum environ- 
ment thus preventing contamination. 

From the standpoint of conventional scanning tunnel- 
ing microscopy (STM) experiments it is also unique in 
both the use of a highly resistive substrate and the ad- 
dition of a bilayer of Ti/Pt electrodes on top of the sub- 
strate. The device geometry is illustrated in Fig. [TJ The 
first barrier separating the STM tip (not shown) from 
the cluster is a vacuum barrier while the second barrier 



FIG. 1: An illustrative drawing of the device geometry. Elec- 
trodes were predeposited on the top of the Si substrate in the 
configuration shown in the diagram and contacted to the sam- 
ple plate. On the back side of the 500 fim thick Si substrate, 
a uniform bilayer of Ti/Pt was deposited for electrical contact 
to the sample plate. 



between the cluster and the substrate is a thin silicon 
oxide barrier. This double tunnel junction arrangement 
is necessary for the resolution of discrete energy levels 
in these nanostructures. The experimental details of the 
cluster fabrication have been published elsewhere. The 
buffer layered assisted growth involves initially deposit- 
ing four monolayers of Xe onto a cold substrate that is 
held at temperatures less than 50K. The pressure of the 
Xe is controlled upon entry into the deposition chamber 
by a capacitance manometer. The condensed Xe is then 
subjected to a precisely controlled exposure of Pb va- 
por flux. The average film thickness of Pb was less than 
0.2 A as measured by a calibrated quartz crystal oscilla- 
tor. Subsequently the substrate is slowly warmed to room 
temperature such that the Xe desorbs and the Pb clus- 
ters softly land onto the substrate. The resulting sample 
which is held on a rotable liquid helium cooled transfer 
rod is moved in situ into the STM chamber through a 
large gate valve and placed via a wobble stick onto the 
STM stage. All STM measurements were made with the 
STM operating at 4.2 K although the electron tempera- 
ture is significantly larger. The STM tip that was used 
in the experiment was made from tungsten wire. The 
tip was characterized prior to making spectroscopy mea- 
surements by demonstrating atomic scale resolution of a 
graphite surface. 

After fabrication, the Pb clusters were characteristized 
using STM. The heights of the clusters were measured 
and found to be mostly between 3 A to 12 A, with an 
average cluster height of 8 A. When a bias voltage was 
applied between the tip and counter-electrodes on top 
of the Si substrate with the cluster in between, the cur- 
rent was measured through the cluster. Current flows 
when the levels of the cluster come into juxtaposition 
with the level(s) at the interface state or another set of 
discrete energy levels from a neighboring cluster. The 
peaks observed in the current-voltage curves (Fig. [2]) are 
indicative that resonant tunneling processes are operative 
whenever one of the cluster's quantized states is probed. 
This will be further elaborated upon in the next section. 
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FIG. 2: Tunneling current versus voltage at T = 4.2 K. Tun- 
neling is from a tunsten STM tip into a Pb cluster. Inset: 
30.0 nm x 30.0 nm image of Pb clusters grown by a buffer 
layer assisted growth technique. This image was obtained us- 
ing a bias voltage of -3.0 V with a tunneling current of 2.0 x 
10-^ A. 



III. CURRENT- VOLTAGE CHARACTERISTICS 

In this section we shall be concerned with measure- 
ments of current-voltage characteristics of a series com- 
bination of two sets of discrete energy levels of systems 
separated by insulating barriers. Electron transmission 
through the clusters was probed by the STM tip in spec- 
troscopy mode. The tunneling current was recorded as a 
function of applied sample bias voltage V while the STM 
feedback loop was open. In this case, the tip was held at 
a fixed height above the cluster while the current-voltage 
characteristics were recorded. The data presented in this 
work were all obtained in voltage step sizes of 480 fiV or 
612 fiV in a time frame > 0.2 ms per point. 

Due to resonant tunneling between the tip, the cluster, 
the interface state, and finally the Ft electrode, the I — V 
characteristics exhibit peaks as shown in Fig. [2l These 
peaks are irregularly spaced and varied in height and are 
dependent on cluster size. The statistics of peak heights 
and spacings will be discussed in the next section (Sec. 
IV). Several hundred I—V curves were obtained on each 
cluster and their character varied along a given cluster in 
a nonsystematic manner. 

The magnitude of the differential conductance is an ap- 
propriate measure of whether or not the tunneling regime 
is the correct description of the transport. This will 
be addressed in this section along with a brief discus- 
sion of elastic and inelastic tunneling as it pertains to 
the linewidths. Also the relationship between the I — V 
curves and the sizes of the clusters will be discussed. But 



first we will take up the issue of resonant tunneling in this 
system. The data is quite different from that which would 
be obtained in a geometry in which metallic atoms (or 
clusters) are deposited on an oxidized metallic surface 
and transmission is via nonresonant processi^ 



A. Resonant Tunneling Mechanism 

The observation of peaks instead of steps strongly sug- 
gests that the transport through the cluster is due to 
resonant tunneling processes.—"— One possible expla- 
nation of how this occurs is the existence of an interface 
state between the Fb clusters and the Si substrate. It 
is known from photoemission experiments that such an 
interface exists between Fb and Si(lll) and this inter- 
face state is nearly dispersionless<^ While there are dif- 
ferences in resistivity between the Si substrate that was 
used in our study and the photoemission experiment, we 
assume that such a state exists. 

In all the I — V spectra, the tunneling current is sup- 
pressed around zero bias followed by a series of peaks 
at both negative and positive bias. The current is sup- 
pressed when the energy levels of the cluster are not in 
registry with the interface state. The double barrier tun- 
nel junction that is realized by positioning the STM tip 
over a Fb cluster is highly asymmetric with the first tun- 
neling barrier between the tip and the cluster and the sec- 
ond barrier between the cluster and the substrate. This 
asymmetry is manifested in the asymmetry of the I—V 
curves about zero bias. 

Resonance occurs when the Fermi wavelength spans 
the length of the cluster. The Fermi wavelength for Fb 
is 4 A at room temperature.-^- Given that the average 
height is 8 A, and measurements were carried out at liq- 
uid helium temperatures, the assumption that the tun- 
neling mechanism is a resonant process is appropriate. 



B. Level Spacing as a function of Cluster Size 

Table U catalogs the Fb clusters with the estimated 
number of atoms per cluster, and estimated and mea- 
sured mean level spacings. The volume of the clusters 
was calculated assuming that each cluster is a hemispher- 
ical cap such that the volume. Vol. is ^(3r^-|-/i^), where 
r is the radius and h is the height of the cluster. Since 
the clusters were irregularly shaped, their radius was es- 
timated from r = \ — where I is the length of the cluster 

y TT ^ 

and w is the width. The estimated number of atoms was 
calculated from Vol./ {\cl^^ where a is the lattice pa- 
rameter of Fb and is equal to 4.95 A^SS, The estimated 
mean level spacing, which was calculated using the nearly 
free electron model, is (A) = J^^^ y^; and compared to 
the measured mean level spacings. Several of the clus- 



ters' measured mean level spacing were in the range of 
8— lOmey. There are several explanations for this. First, 
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the majority of the clusters had the same height which 
is a result of the buffer layer assisted growth technique 
when using a four monolayer thick buffer layer. Second, 
the calculation of the volume is a rough estimate and not 
entirely accurate in that it oversimplifies the actual shape 
of the clusters. Third, the clusters were not well isolated 
from each other thus the width and the length could be 
larger if more than one cluster was taken into consider- 
ation. Also, tip convolution effects obscure the actual 
width and length of the clusters. Nonetheless there is an 
apparent dependency on the level spacing energy on the 
size of the clusters. 



C. Lineshapes and Differential Conductance 

Information about the intrinsic lifetime of the elec- 
tronic states of a single metal cluster and the tunneling 
process (inelastic or elastic) is in part difficult to obtain 
because of the temperature of the reserviors. This can 
be circumvented by having two sets of electronic states 
separated by a barrier through which the current is mea- 
sured. In this tunneling scenario, electron transport is 
expected to occur only when the energy levels of the 
two states are aligned with one anothei— and the up- 
per bounds of the lifetimes of energy levels of varying 
size clusters can be estimated from the average widths of 
the resonances that results from these level alignments. 
The reason that these widths are upper bounds is that 
they depend on both the reciprocals of the lifetimes of 
the electronic states and the tunneling rates. We fit the 
full width half maximum of the peak widths in the I- 
V curves for three clusters to Lorentzians (Lorentzians 
provided the best fit to the data) and found that they 
varied depending on the cluster size. For cluster 1, the 
average peak width was 3.54 meV (tunneling time: 1.9 
X 10"^'^ seconds), for cluster 7, the average peak width 
was 3.06 meV (tunneling time: 2.4 xlO""'^^ seconds) and 
for cluster 12, the average peak width was 1.06 meV (tun- 
neling time: 6.2 xlO~^^ seconds) indicating that it takes 
longer to tunnel through a larger cluster than a smaller 
one. Lorentzian peak widths are also an indication that 
the tunneling process is elastic. Moreover, the absence 
of diffusion (the transport is ballistic) together with the 
spatial resolution of the STM enables detailed studies of 
the clusters' electronic states. 

Figure [3] shows a plot of the differential conductance in 
units of e^/h plotted against positive bias for one of the 
clusters. This is supporting evidence that the features are 
related to elastic tunneling processes as the conductance 
is less than one and the peaks are Lorentzians. 



D. Position Space Representation of Eigenstates 

An important feature of the measurements was that 
the current-voltage characteristics varied with position 
on a cluster as shown in Fig. [D The peak heights. 
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FIG. 3; Differential conductance (in terms of e^/h) versus 
voltage at T= 4.2K for cluster no. 12. The symmetry about 
the y-axis and the magnitude of the peaks indicates that elas- 
tic tunneling processes are relevant. 



the number of peaks, and voltages where the peaks oc- 
curied varied randomly as a function of position. This 
dependence is reminiscent of a quantum chaotic system 
in which the amplitude of the eigenfunctions are known 
to vary in a complicated patterni^l There was no notice- 
able systematic variation in the curves as a function of 
distance. It was not possible to reproduce the curves be- 
cause of the difficulty in repeating measurements at the 
exact same pixel point due to drift in the piezos which 
control the tip's position. 

For a given current profile across a cluster, a map of the 
current as a function of the voltage bias can be obtained 
as shown in Fig. [5] For cluster no. 12, such a map dis- 
plays strong clustering of current peaks as a function of 
applied bias voltage although there is no obvious voltage 
dependence. 



E. Other possible tunneling scenarios 

The current voltage characteristics bare a strikingly 
similar resemblance to the current vs gate voltage curves 
seen in the results reported by Kouwenhoven, Austing 
and Tarucha.^® It might be the case, that neighboring 
clusters are actually gating the cluster under investiga- 
tion. This would offer an alternative explanation to the 
position dependence of the I-V curves as the capacitance 
will change due to the distance between the gating clus- 
ter and the site on the cluster that was being probed by 
the STM tip. Further work is needed to elucidate the 
mechanism behind the peaks in the current-voltage char- 
acteristics. One such possible experiment would be to 
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Cluster 


Length 


Vol 


<N> 


Est. <A> 


Meas. <A> 


dE drso 


Index No. 


(nm) 


{nmf 




(meV) 


(meV) 


(ft= 1) 


1 


2.7 


2.6 


85 


36.8 


24.2 ± 5.39 


.083 


2 


3.32 


3.81 


125 


25.1 


6.35 ± 0.74 


.070 


3 


2.29 


3.86 


127 


24.76 


9.02 ± 1.45 


.047 


4 


3.61 


6.29 


207 


15.2 


10.37 ± 0.43 


.046 


5 


3.13 


6.3 


208 


15.15 


9.41 ± 1.12 


.04 


6 


3.04 


7.36 


243 


12.99 


9.13 ± 0.29 


.033 


7 


2.7 


7.4 


244 


12.91 


10.8 ± 0.54 


.029 


8 


3.88 


7.8 


257 


12.25 


8.39 ± 0.28 


.04 


9 


5.0 


8.2 


270 


11.65 


9.01 ± 0.26 


.049 


10 


3.73 


9.62 


317 


9.93 


8.58 ± 0.34 


.031 


11 


3.12 


10 


329 


9.55 


5.37 ± 1.53 


.025 


12 


4.3 


12.7 


418 


7.52 


9.6 ± 1.25 


.027 


13 


2.54 


13.2 


435 


7.23 


8.98 ± 1.63 


.0154 


14 


5.49 


19.4 


638 


4.94 


7.35 ± 1.30 


.023 


15 


4.12 


21.5 


709 


4.44 


7.49 ± 0.37 


.015 


16 


4.76 


27.1 


892 


3.53 


4.73 ± .87 


.014 



TABLE I: Summary of physical parameters of Pb clusters of different shapes and sizes. As discussed in the text, the Vol. is 
the volume of a cluster and is estimated assuming it is a hemispherical cap. The number of atoms, jN/,, and the est. mean 
level spacing are calculated from this est. volume. The est. mean level spacing is based on the nearly free electron model and 
compared to experimentally measured mean level spacing values. The expression dE drso is related to the cluster's size and is 
described in the text (see Sec. IV D). 



fabricate the clusters without the electrodes configured 
on top of the substrate. This would provide a more uni- 
form distribution of clusters and increase the separation 
between clusters since the clusters tended to conglomer- 
ate near the electrodes. 



IV. STATISTICAL DISTRIBUTIONS 

Current peaks in the scanning tunneling spectroscopy 
were not expected. However, in an attempt to interpret 
the results the data was analyzed with the assumption 
that the peaks were signatures of discrete energy levels. 
Exact agreement with different sets of discrete energy lev- 
els for both the negative and positive bias voltages would 
demand that the capacitances between both the tip and 
cluster and cluster and substrate are substantially dif- 
ferent, which we believe is the case. Regardless, from 
the I-V measurements themselves, it is not possible to 
confirm that the peaks in the negative and bias voltages 
arise from different energy levels without further analy- 
sis. Therefore, based on the premise that the discrete 
energy spectra in metallic clusters should follow RMT, 
we studied the statistical distributions of both the peak 
spacings and peak heights. Peaks in the I-V characteris- 
tics were identified after running a smoothing program^S 
twice through the data. This program was used to iden- 
tify all the peaks. Since the voltage range over which 
peaks were found in this study was small (-35 meV to 
-1-30 meV) , the likelihood of missing any peaks was elim- 



inated. A practical matter connected with the reliability 
of the histograms must also be mentioned. The number 
of bins used in the histograms was the range divided by 
the experimental step size. The histograms were rescaled 
using the measured mean level spacing to normalize the 
mean to unity. The fits to the resulting histogram were 
carried out using MINOS^S which is a minimization algo- 
rithm implemented in MINUITi^ The parameters that 
are obtained are those that correspond to a minimum chi- 
squared value. Also, the peak widths are approximately 
ten times greater than the bin size, thus the peaks spacing 
distribution that was generated is not due to a " binning" 
effect. 



A. Wigner-Dyson Statistics 

The histogram of normalized peak spacings (normal- 
ized to the mean spacing of each individual trace) was fit 
by the distribution function 

P(s) = bfiS^ eyii>{-ap s"") (1) 

The distribution function P(s) is a probability density 
defined such that the area under the curve is one. Here 
the normalized mean spacing, s, is simply A / (A) , with 
A representing the level spacing and (A) the mean spac- 
ing. Equation [T] can represent the orthogonal (/? = 1), 
unitary {(3 = 2) and symplectic (/3 = 4) ensembles^ 
that correspond to processes with different symmetries. 



5 




FIG. 4: STM image (12 nm x 8.0 nm) of Pb clusters. The 
I — V curves correspond to three different positions along a 
cluster. There is not any apparent correlation between the 
curves. The current setpoint = 1.0 nA. 



The orthogonal case corresponds to time reversal sym- 
metry being preserved in the absence of a magnetic field 
and describes the results presented here. In the statisti- 
cal analysis of this histogram fits by Wigner-Dyson, Pois- 
son, Gaussian, and Lorentzian distributions were made. 
From the values of x^(iiot shown) it is clear that the 
Wigner-Dyson distribution provides the best fit to the 
data with — n/A, bp — 7r/2 and cp = 2. In Fig. [HI the 
histogram of peak spacings for this cluster, showing the 
Wigner-Dyson and Poisson fits is plotted. (In the figure, 
the fitted Poisson distribution is not one about the ori- 
gin since the parameters were allowed to float in order to 
minimize chi"2). 



B. Porter-Thomas Statistics 

The following form^^. 



P{I) 



<I> 



exp 



< / > 



(2) 



was fit to the data, where / is the peak current and 
< I > the mean peak current. In this analysis, param- 
eters specific to the Porter-Thomas and Poisson distri- 
butions, which were deemed relevant, were used. The 
Porter-Thomas distribution (with a = {2tt)^2^ b = —1/2 
and c — 1/2) provided a somewhat better fit to the data 
than the Poisson distribution. Figure [7] shows a plot of 
the histogram along with curves associated with the best 
fits of the Porter-Thomas and Poisson distributions. The 
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FIG. 5: A map of the current as a function of position along 
a chaotic cluster (catalog number 4) as the bias voltage is 
varied. There is noticable clustering of the current peaks. 
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FIG. 6: Histogram of peak spacings for cluster no. 9. The 
solid curve is the fit for the Wigner-Dyson distribution. The 
dotted line represents the fit for the Poisson distribution. 
There are 413 peak spacings that comprise this histogram 
normalized to the mean voltage spacing. 
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FIG. 7; Histogram of peak heights for cluster no. 9. The solid 
curve is the fit for the Porter- Thomas distribution. The dot- 
ted line represents the fit for the Poisson distribution. There 
are 851 peaks heights that comprise the histogram. 



results of this analysis support the interpretation that 
these measurements are yielding spectroscopic informa- 
tion relating to the energy levels. 



while the bottom histogram is strongly Wigner-Dyson 
like. This is in agreement with other experimental situa- 
tions, such as the He atom in a magnetic fieldi^ and the 
acoustic resonances of Al blocks where deformations in 
the shape of the block away from parallelpiped generate 
partially chaotic acoustic wavesj^ Completely chaotic 
systems are rare and difficult to achieve experimentally. 

The modified Berry-Robnik equation^ indicates what 
percentage of the distribution is Poisson like and what 
percentage is Wigner-Dyson like and accounts for the 
small level repulsion that is observed in most systems. (It 
is a modified version of the Berry-Robnik distribution— 
which works well for mixed chaotic and regular systems 
at the tail of the distribution, but fails at small level 
spacings.). Thus the modified Berry-Robnik equation is 



(3a) 



P(s) cx q^F(^)e-«''erfc (1 - g) s 



|(l-g)^ + 2g^^(^)](l-g) 



(3b) 



where erfc is the complimentary of the error function and 
is defined as 



erfc {x) =l-erf (x) — —= e " du 



(4) 



In addition, F{x) is defined as follows: 



F{x) = l- 



(5) 



When q = the Wigner-Dyson term is retrieved and 
when q=l, the distribution is a Poisson distribution. 
Equation [3] describes a mixed state between these two 
limits and addresses the crossover regime in which the dy- 
namics are partially chaotic. In Table [Til the fitted values 
for the modified Berry Robnik equations are summarized 
for different sized Pb clusters. 



C. Crossover between Poisson-like and Wigner- 
Dyson Statistics 

Experimentally one should observe variations in the 
distributions of the energy levels based on the cluster's 
shape. Keeping all the experimental parameters the 
same, i.e. the same tip-cluster height distance and the 
same voltage scale, the distributions of the eigenvalues 
change depending on the cluster's shape. In the work pre- 
sented here, the majority of distributions for the different 
clusters fell into the regime between Wigner-Dyson and 
Poisson statistics, which describes a partially chaotic sys- 
tem. This is plotted in Fig. [S] where clusters of the same 
approximate volume but have different shapes are consid- 
ered. The top histogram in Fig. [5] is more Poisson-like 



D. Absence of Charging Energy, Superconductiv- 
ity and Spin-Orbit Scattering 

An important issue is the role of charging energy in 
the proposed two-step tunneling process. It is likely that 
the charging energy is absent because of additional ca- 
pacitance arising from the proximity of neighboring clus- 
ters. This will lower the charging energy significantly so 
that the peaks that are observed are only those related 
to the eigenstates of the clusters. Absence of charging 
energy was also recently reported from x-ray photoemis- 
sion spectroscopy measurements— on gold nanocrystals 
which were self-assembled using wet chemistry. In this 
study, as the mutual separation between gold clusters 
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FIG. 8: Transition to chaos in the level statistics of Pb clus- 
ters of the same approximate volume (« 10 nm^) but with 
different shapes. The solid black curve is the Wigner-Dyson 
distribution, while the dotted black curve is the Poisson dis- 
tribution for comparison. The cluster's catalog number is in 
the upper right hand corner. 



Cluster 


#1-V 
curves 


# Peak 
Spacings 


q 


V 


1 


<10 


<30 


— 





2 


180 


35 


— 





3 


119 


137 


.59±0.71 


.65 ±0.22 


4 


137 


144 


0.02 ±0.06 


— 


5 


125 


16 


— 


— 


6 


158 


254 


-3.1e-9 ±4.39e -02 


— 


7 


<10 


<30 


— 


— 


8 


207 


273 


.14±0.085 


.45±0.17 


9 


280 


413 


0.02±0.037 


— 


10 


131 


181 


.06 ±.062 


.18±0.25 


11 


316 


142 


1.37 ±0.56 


.73±0.06 


12 


<10 


<30 


— 


— 


13 


193 


237 


0.004±0.05 




14 


405 


64 


1.28±0.58 


.7±0.16 


15 


220 


105 


-0.0086±.07 




16 


263 


450 


0.25±0.07 


.48±.08 



TABLE 11: Summary of the statistics of different sized Pb 
clusters. When q=0 or close to it, WD statistics prevails and 
the value of v becomes meaningless, and when q = 1 Poisson 
statistics is the appropriate description of the distribution. 
(The low statistics for cluster numbers 1, 2, 5, 7, and 12 
resulted in undefined q and v values.) 



was reduced to distances less than Inm by varying the 
Ugand size, the charging energy was eliminated. New ex- 
perimental conditions need to be reahzed in our work in 
which the clusters are well separated from each other in 
order to resolve this issue about the charging energy. 

The apparent absence of superconductivity is a result 
of the mean level spacing exceeding the energy gap of 
superconducting Pb which is ~ 2 meV . For this reason 
one would not expect to observe features in the tunneling 
characteristic associated with superconductivity.'^- This 
was also observed in the pioneering work by Ralph, Black 
and Tinkham.'^ 

Lead is a strong spin-orbit coupling material and the 
appropriate description to describe strong spin orbit cou- 
pling is a Gaussian Sympletic Ensemble (GSE) instead of 
Gaussian Orthogonal Ensemble (GOE or Wigner-Dyson) 
distribution. How strongly spin orbit coupling adheres 
to GSE instead of GOE can be related to the cluster's 
size^ as follows: If Wigner-Dyson (GOE) statistics pre- 
vail, then the following criteria based on the uncertainity 
principle (h = 1) should be met: 6E6tso ^ 1- Like- 
wise, if the appropriate statistics is related to Gaussian 
Sympletic Ensemble (GSE), then the following criteria is 
applicable: SEStso <C 1 . In these relations, SE is the 
mean energy level spacing, or one divided by the density 
of states. The term Stso is the time that it takes to 
flip a spin or 6tso = ^ where L is the length of the 
cluster, and w/ is the Fermi velocity. This assumes that 
the lateral transport in the dot is ballistic with the only 
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FIG. 9; Current as a function of position on a chaotic cluster 
(cluster no. 4) at 26 meV. 



scattering occuring at the boundaries. Thus the expres- 
sion SEStso is equal to p^?^. In the Table |T] the values 

for SEdrgo are calculated based on cluster size instead of 
the actual mean level spacing which is smaller than the- 
oretical estimates. This upper limit indicates that the 
appropriate distribution should be GSE instead of GOE. 
This contradicts our experimental findings. 

It should be noted that in studies of Au grains (<5 
nm), it was found that the spin-orbit scattering was 
suppressed.-^ It was speculated that the origin of the sup- 
pression of the spin-orbit scattering was the granularity of 
the weakly coupled grains. Estimates of this suppression 
were carried out by making weak localization measure- 
ments. This involves a study of the orbital effect in an 
applied magnetic field. While we could not make such 
estimates (our apparatus lacks a magnetic field) using 
similar arguments, we may have observed a suppression 
of spin-orbit scattering (based on the statistics). 

V. IMAGING QUANTUM CHAOS WITH AN 
STM 

Finally we turn our attention to imaging the eigen- 
functions of the clusters using a scanning tunneling mi- 
croscope. The wavefunction moduli are directly related 
to the tunneling current, so by mapping the tunneling 
current amplitude at a fixed bias, it is possible to study 



the significance of chaos as a function of cluster shape. 
An image of the current peaks is shown in Fig. [5] for 
cluster no. 4. The current peaks change as a function 
of applied bias (not shown) . 

The current dependence as a function of bias voltage 
can be used to "image" the chaotic nature of irregu- 
larly shaped metallic clusters. This gained prominence 
in the experimental community with the pioneering ex- 
periments of Sridhar and co-workers4i Besides acquiring 
images of the clusters, it is possible to produce images 
of chaotic behavior in the same manner than was carried 
out by Sridhar using electric field to probe the distri- 
butions of eigenfunctions in a "Sinai and rectangular" 
structure. In essence scanning tunneling spectroscopy is 
a technique that can be used to map out the amplitude of 
the eigenfunctions which are otherwise difficult to access. 

VI. CONCLUDING REMARKS AND OPEN 
QUESTIONS 

The empirical evidence indicates that discrete energy 
levels are being accessed from a resonant tunneling ge- 
ometry and that the distributions of these energy levels 
reflect the underlying classical dynamical nature of an 
electron under confinement of its boundary. This means 
that our study of the distribution of spacings gives in- 
formation about the symmetry properties of the physi- 
cal system in question. In this study of metal clusters, 
there are remaining questions that still need be resolved. 
First what is the nature of the tunneling process? Is the 
electron's route through an Si interface state, or through 
neighboring clusters? Future experimental work needs to 
be carried out in which the clusters are separated from 
one another in order to resolve this issue. The second 
question that needs to be answered is why are the spin- 
orbit interactions suppressed? That is, why not Gaus- 
sian Sympletic Ensemble statistics instead of Wigner- 
Dyson statistics? This issue needs to be explored both 
theoretically and experimentally in greater detail. 
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